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ABSTRACT 


>/•  i 


An  approach  was  developed  by  the  author  to  estimateQhe 
effect  of  seakeeping  on  ‘^^h^mission  effectiveness  of  a 
feasible  design^-  The  key  to  this  approach>is  'Uie  applica¬ 
tion  of  a  model  of  the  mission  as  a  system  whei^n  the  func¬ 
tions  which  perform  the  mission  have  the  same  relationship 
to  each  other  as  the  components  do. in  a  physical  system. 

Th^ ^egi^j^ation  in  ‘^h^  performance  of  the  mission  is  calcu¬ 
lated  in  the  same  manner  as ^Che^rel lability  of  a  system  made 
up  of  separate  components.  This  model  is  incorporated  into 
the  feasibility  study  stage  of  the  ship  design  process,  and 
the  modified  procedure  is  described. 


The  approach  is  illustrated  by  applying  it  to  a  single/"’''  - 
example, a  monohul  1  •. ■  'Fhe;^'  f-unctions  necessary  to  perform  a 
given  mission  and  'th«>>relatiOTships  between  separate  func¬ 
tions  were  determined.  Physical  subsystems  which  perfcsrm 
those  functions  were  selected  and  substituted  into  the  func¬ 
tional  relationship.  A  regressi on— based  computer  ship  syn¬ 
thesis  program  for  destroyer— type  ships  was  used  to  obtain 
tb^  hull  form  characteristics  necessary  to  create  a  proto¬ 
type  hull  form.  A  prototype  hull  was  developed  using  a  hull 
form  generation  program  which  uses  parametric  cubic  splines 
to  build  a  definition  of  the  ship's  hul  1 Fi ve  locations  on 
the  hull  were  selected  and  each  subsystem  assigned  to  be  at 
or  near  one  of  those  locations.  The  vertical  motions  of 
those  locations  were  calculated  for  the  case  of  long  crested 
head  seas  using  Lewis-form  approximation  to  the  ship's  hull. 
The  motions  were  used  to  determine  the  relative  effective— 
eness  of  each  subsystem  from  its  performance  degradation 
function.  The  relative  mission  effectiveness  was  computed 
using  the  alg^ithm  used  to  compute  the  reliability  for  a 
subsystem  composed  of  different  components. 
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Title:  Professor  of  Ocean  Engineering 
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CHAPTER  1 


BACK8R0UND 

DEFINITION  OF  SEAKEEPING 

Increasingly,  the  United  States  Navy  is  including 
seakeeping  as  a  per-f ormance  requirement  in  the  design  of 
neM  US  Navy  ships.  The  definition  of  seakeeping  widely 
accepted  by  the  US  Navy  design  community,  and  which  will 
be  used  here  was  given  by  Vice  Admiral  R.  E.  Adamson: 

"Seakeeping  as  it  pertains  to  the  U.S.  NAVY,  is 
the  ability  of  our  ships  to  go  to  sea  and  success¬ 
fully  and  safely  execute  their  mission  despite 
adverse  environmental  factors."  HI 

The  ship's  response  to  waves  causes  a  decrease  in  the 
ability  of  equipment  and  personnel  to  perform  many 
functions.  Seakeeping  performance,  also  called  mission 
effectiveness,  is  the  measure  of  the  ship's  ability  to 
carry  out  a  mission  fully  and  completely.  More 
specifically,  a  ship  with  ideal  seakeeping  qualities  is 
capable  of  performing  its  mission  in  all  but  the 
roughest  seas;  it  is  not  limited  to  a  particular 
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direction  with  respect  to  waves,  nor  to  less  than  its 
rated  top  speed,  nor  is  any  operation  it  may  be 
Assigned  to  conduct  limited  by  sea  state  alone. 

CURRENT  PRACTICE  IN  DESIGN  FOR  SEAKEEPING 

The  Naval  Sea  Systems  Command  (NAVSEA)  is  the  chie-f 
design  agent  -for  new  ships  -for  the  U.S.  Navy.  At 
NAVSEA,  efforts  are  focused  on  three  design  areas  when 
designing  for  good  seakeeping  qualities.  123  The  first 
is  the  hull  configuration.  Within  the  limits  of 
sever -:«1  non-seakeeping  constraints  which  deal  with 
power  and  speed,  arrangements,  center  of  gravity, 
stability  and  structural  strength,  the  designer  strives 
for  the  shape  which  minimizes  overall  ship  motions  and 
events  such  as  slamming  and  deck  wetness.  The 
designer  tries  to  optimize  features  such  as  freeboard, 
shear,  and  secondary  hull  form  parameters  133.  The 
second  area  is  the  arrangement  of  motion  sensitive 
equipment  in  order  to  minimize  the  effects  of  motion. 
The  third  area  is  the  design  of  the  equipment  to 
minimize  the  effects  of  the  ship’s  motion  on  the 
performance  of  the  equipment.  In  the  third  area 
equipment  is  usually  not  designed  with  a  particular 
ship  in  mind,  but  with  the  intention  of  installing  it 
on  many  classes  of  ships.  In  the  first  two  areas  the 
optimization  is  done  once  a  design  alternative  has  been 
selected  for  further  development. 


9 


AN  ALTERNATIVE  APPROACH 

The  optimization  of  the  design  for  good  seakeeping 
qualities  begins  after  the  feasible  design  process  has 
selected  "best”  alternative.  There  is  currently  no 
method  for  predicting  quantitatively  how  the  seakeeping 
qualities  of  a  design  alternative  will  affect  its 
mission  effectiveness.  If  an  estimate  of  the 
degredation  of  the  mission  effectiveness  of  a  ship  due 
to  seakeeping  could  be  make  for  alternatives  during  the 
feasible  design  stage,  this  would  give  the  designer  the 
ability  to  include,  even  though  to  a  very  limited 
extent,  the  effects  of  seakeeping  in  the  selection  of  a 
"best"  alternative  design. 

This  thesis  presents  one  approach  to  estimating  the 
mission  performance  degr^ation,  termed  mission 
effectiveness,  and  proposes  a  modified  feasible  design 
process  incorporating  this  approach.  The  author  shows 
how  the  modified  process  could  be  applied  to  a 
monohull,  and  illustrates  the  application  with  a 


numerical  example 


CHAPTER  2 


A  DEFINITION  OF  MISSION  EFFECTIVENESS 


A  measure  o-f  a  ship's  seakeeping  ability,  as 
seakeeping  is  defined  above,  is  a  measure  of  how 
effectively  a  mission  is  performed  on  the  open  ocean  in 
the  presence  of  waves  and  other  natural  environmental 
factors.  Relative  mission  effectiveness  is  the  ratio  of 
the  ship's  ability  to  carry  out  a  mission  in  a 
specified  sea  state  to  the  ships  ability  to  carry  out 
that  same  mission  in  calm  water.  This  will  generally 
be  referred  to  as  mission  effectiveness  or  performance. 
The  specific  criteria  for  measuring  effectiveness  will 
vary  from  mission  to  mission.  In  each  case,  though, 
the  relative  effectiveness  is  a  dimensionless  ratio. 

A  PROPOSED  MODEL  FOR  A  MISSION  STRUCTURE 

A  mission  is  a  complex  operation  which  can  be 
described  in  simplest  terms  by  a  single  word  or  short 
phrase  such  as  "transport  cargo",  "transport 
passengers",  "conduct  anti-submarine  warfare",  or 
"conduct  search  and  rescue".  These  complicated 
operations  can  be  broken  down  into  several  simpler 
operations,  or  functions.  Many  functions  may 
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themselves  be  complicated  operations  and  can  be  further 
broken  donnn  into  sub-functions.  Several  functional 
levels  may  be  described  until,  at  some  point  functions 
can  be  related  to  equipment  or  subsystems  which  perform 
them. 

If  a  mission  is  viewed  as  a  system  of  functions, 
and  if  the  functions  are  treated  as  independently 
degraded  by  ship  motions,  then  the  method  used  for 
computing  system  reliability  for  RMA  analysis  can  be 
used  to  calculate  mission  effectiveness  143.  Mission 
effectiveness  is  the  product  of  the  effectiveness  of 
the  different  functions  which  make  up  the  mission.  The 
functions  which  perform  a  mission  are  performed  either 
in  series  or  in  parallel.  For  the  purposes  of 
determining  mission  effectiveness,  a  mission  can  be 
diagrammed  in  the  same  way  a  system  of  components  would 
be  diagrammed  for  the  purpose  of  making  a  system 
reliability  calculation. 

HOW  FUNCTIONS  ARE  RELATED  TO  EACH  OTHER  IN  THIS  MODEL 

Functions  are  performed  either  in  series  or  in 
parallel.  The  actual  determination  of  which  functions 
are  in  series  with  each  other  and  which  are  in  parallel 
are  very  much  a  matter  of  experience  and  good 
judgement.  One  criterion  which  can  help  is  to  ask  if  a 
mission  could  be  performed  if  a  particular  function 
could  not  be  performed  at  all.  If  the  mission  would 


fail,  then  that  function  is  in  series  Mith  other 
functions.  If  the  Mission  could  still  be  carried  out, 
then  the  function  is  in  parallel  Mith  other  functions. 

Several  functions  or  sub~f unctions  say  be  grouped 
together  in  parallel  and  this  group  in  turn  may  be  part 
of  a  serial  group,  or  vice  versa.  Groups  may  be  nested 
inside  of  groups.  Then  in  order  to  make  the  effec¬ 
tiveness  calculations,  the  value  for  the  innermost 
group  must  be  calculated  first.  The  hypothetical 
mission  graphed  in  figure  (2-1)  is  an  example.  The 


FIGURE  (2-1) 


A  HYPOTHETICAL  MISSION  STRUCTURE 


mission  is  represented  by  ten  functions  or  sub- 
functions,  fl  through  fiO.  There  are  five  subgroups; 

PI  consisting  of  functions  f2  and  f3,  P3  consisting  of 
subgroups  SI,  S2,  and  function  f8,  and  SI  consisting  of 
fS  and  P2,  and  so  on. 

SUBSYSTEMS  AND  THEIR  RELATION  TO  FUNCTIONS 

Actual  performance  of  fiuictions  is  done  by  sub¬ 
systems,  equipment  and  personnel.  Hereafter,  the  term 
subsystem  will  include  personnel.  Just  as  several 
functions  may  perform  a  mission,  several  subsystems  may 
perform  a  particular  function.  The  same  process  of 
determining  the  series-parallel  relationships  of  the 
functions  in  a  mission  to  each  other  may  be  applied  to 
several  subsystems  which  perform  a  single  function. 

If  the  performance  of  a  subsystem  or  piece  of 
equipment  is  degraded  by  the  motions  of  the  ship,  the 
effectiveness  of  the  function  performed  by  that 
particular  subsystem  is  also  degraded.  The  mission 
effectiveness  is  also  degraded.  The  definition  of 
relative  mission  effectiveness  introduced  at  the 
beginning  of  this  chapter  can  also  be  applied  to  the 
performance  of  subsystems.  The  relative  performance  of 
a  subsystem  is  the  ratio  of  its  actual  performance  when 
degraded  by  ship  motions  to  its  performance  under  calm 
water  conditions. 

A  performance  degredation  function  is  the 
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correlation  of  this  relative  performance,  or  degraded 
performance  to  some  significant  measure  of  the  motion 
which  causes  it.  The  information  that  is  available 
currently  for  the  performance  degredation  of  existing 
subsystems  is  sketchy  at  best.  It  is  often  based  on 
the  only  observations  available,  the  reports  and 
comments  of  the  operators.  About  the  only  effectiveness 
carefully  studied  is  human  performance,  which  is 
significantly  nonlinear.  Other  studies  have  focused  on 
complex  functions,  such  as  fuelling  at  sea  which  has 
many  steps  in  its  execution.  Many  other  performance 
degredation  “functions"  consist  of  two  points  on  a  plot 
of  percent  effectiveness  versus  some  motion  amplitude. 
One  point  is  the  largest  reported  amplitude  with 
"satisfactory”  or  fully  effective  performance,  the 
other  is  a  reported  amplitude  where  the  performance 
could  not  be  performed  at  all.  A  straight  line  between 
these  two  points  constitutes  the  best  guess  (and  this 
is  recognized  and  openly  admitted  as  a  best  guess)  C51 
as  to  the  performance  degredation  of  a  particular 
activity  or  function.  Efforts  are  underway  to  gather 
more  accurate  data  and  to  formulate  better  ways  to 
gather  data  C51. 

CALCULATION  OF  RELATIVE  MISSION  EFFECTIVENESS 

Calculation  is  performed  by  taking  the  product  of 
the  quantitative  values  of  a  function’s  effectiveness. 


15 


For  functions  in  series: 
n 

E  =  YT 
i*l 

and  for  functions  in  parallel: 
n 

E  =  1  -  I  (1-ei)  (2-2) 

i  =  l 

Once  the  values  of  effectiveness  for  each  individual 
function  has  been  determined,  the  value  of  the 
effectiveness  of  any  subgroups  must  be  determined 
before  the  effectiveness  of  the  mission  can  be 
calculated. 

Referring  to  the  hypothetical  mission  in  figure 
(2-1)  as  an  example,  first  the  effectiveness  for 
functions  fl  through  flO  must  be  calculated  separately 
Then  the  effectiveness  for  blocks  P2  and  S2  are 
calculated.  Let  Vj  be  the  value  of  individual 
effect! venesses : 

Vp2  =  1  -  (l-Vf^) (l-Vfy)  and 
Vs2  =  V^9''fl0 

Next  the  effectiveness  for  block  SI  can  be  computed: 

Vsi  =  Vf5Vp2 

Now  the  effectiveness  for  the  block  P3  can  be  computed 
Vp3  =  1  -  (l-Vsi) (l-VfQ) (1-Vs2) 


Vpi  is  likewise  calculated,  then  the  relative  mission 


CHAPTER  3 


A  MODIFIED  FEASIBLE  DESIGN  PROCESS 


The  overall  design  process  is  composed  o-f  several 
phases  o-f  which  the  -feasible  design  stage  is  the  -first. 
The  goals  at  this  stage  oT  the  design  process  are  to 
-first,  eliminate  designs  and  concepts  which  will  not 
meet  the  requirements  and  constraints,  and  second, 
choose  that  alternative  which  will  best  meet  the 
requirements  and  still  remain  within  the  limits  set  by 
the  constraints. 

The  -feasible  design  process  proposed  here  is  shown 
graphically  in  -figure  (3-1).  A  brief  synopsis  follows. 
Mission  requirements  are  considered  a  given  input.  The 
necessary  functions  which  support  and  carry  out  the 
mission  must  be  derived  from  the  requirements.  The 
series  -  parallel  structure  of  the  functional 
relationships  must  be  developed.  Equipment  and 
subsystems  must  be  selected  from  what  is  available  or 
expected  to  be  available  when  the  ship  is  produced. 
These  subsystems,  along  with  performance  criteria  form 
the  input  to  a  ship  synthesis  model  which  takes  the 
payload  and  performance  criteria  and  calculates  a 
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FIGURE  (3-1) 

PROPOSED  FEASIBLE  DESIGN  PROCESS 
This  representation  emphasizes  where  the  estimation  o-f 
seakeeping  performance  fits.  It  does  not  show  all  details 
of  the  design  process. 
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preliminary  estimate  of  the  ship's  physical 
characteristics.  A  hull  form  for  the  ship  is  developed 
and  its  motion  response  to  waves  calculated.  The 
performance  degredation  of  the  subsystems  resulting 
from  the  ship's  motions  is  found  and  the  mission 
effectiveness  is  estimated.  Alternative  hull  types  and 
subsystems  will  give  rise  to  several  alternatives  which 
are  evaluated  in  the  same  way.  Their  performance, 
along  with  how  well  the  alternatives  meet  other 
non— seakeeping  criterion  are  compared,  and  the  one  that 
appears  best  is  selected  for  further  design 
devel opment . 

MISSION  REQUIREMENTS 

The  process  for  designing  a  new  ship  in  the  U.S. 
Navy  begins  with  setting  requirements  for  force  levels 
based  on  national  foreign  and  strategic  policy.  The 
mission  requirements  for  a  ship  are  then  based  upon 
where  this  ship  design  is  expected  to  fit  in  our  naval 
forces.  The  Mission  Requirements  must  next  be  analyzed 
and  mission  areas  defined.  Specific  mission  areas  may 
be  listed  as  part  of  the  mission  requirements. 
Alternatively,  the  mission  requirements  may  indicate 
that  the  new  design  is  to  replace  or  be  similar  to  an 
existing  class  of  ships. 

DETERMINATION  OF  FUNCTIONS  TO  PERFORM  THE  MISSION 

The  first  step  in  determining  the  effectiveness  of 
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a  mission  is  to  determine  what  separate  steps,  or 
-functions  are  performed  in  carrying  out  the  mission.  A 
standard  list  of  mission  areas  which  must  performed  by 
U.S.  Navy  ships  is  a  part  of  OPNAV  3501. 2E(U), 
reference  C63.  This  also  subdivides  the  mission  areas 
into  the  capabilities  necessary  to  carry  out  these 
•nissions.  In  light  of  the  extensive  breakdown  of 
missions  into  their  necessary  functions  which  it 
provides,  the  determination  of  functions  should  be 
straightforward.  If,  for  some  reason  determining  the 
necessary  functions  is  not,  functional  flow  diagramming 
techniques  Ca]C91  can  also  be  used  to  identify  the 
individual  functions  necessary  to  carry  out  the 
mission,  and  the  relationships  between  functions.  A 
top  level,  or  gross  operational  level  functional  flow 
diagram  is  developed  from  the  mission  requirements. 

Then  each  function  can  be  broken  up  into  sub-functions 
and  diagrammed  in  greater  detail,  if  necessary. 
DETERMINATION  OF  THE  MISSION  STRUCTURE 

Once  the  functions  which  perform  the  mission  have 
been  identified,  the  series  —  parallel  relationships 
can  be  developed  as  described  in  chapter  2.  A  diagram 
illustrating  the  functional  relationships  for  each 
mission  may  be  found  to  be  helpful  for  this  step. 
SELECTION  OF  SUBSYSTEMS 

Having  listed  and  diagrammed  the  functions,  the 


actual  subsystBMS  must  b«  s«l acted.  There  may  be 
several  candidate  subsystees  capable  o-f  performing  any 
particular  function.  Several  lists  of  subsystem  which 
make  up  the  various  feasible  alternatives  should  be 
made.  To  illustrate,  suppose  a  mission  has  four 
functions  FI  through  F4.  It  is  determined  that  the 
mission  could  not  be  accomplished  at  all  if  FI  cannot 
be  accomplished,  but  F2,  F3,  or  F4  could  fail 
completely  and  the  mission  could  still  be  performed. 
Furthermore  it  is  determined  that  F4  cannot  be 
accomplished  if  F3  cannot  be.  Thus  F3  is  in  series 
with  F4,  this  pair  is  in  parallel  with  F2,  and  FI  is  in 
series  with  these  as  shown  in  figure  (2-2). 


FIGURE  (2-2) 

ILLUSTRATION  OF  SERIES-PARALLEL  MISSION  STRUCTURE 

Further,  suppose  there  are  two  major  alternative 
hull  types.  In  selecting  subsystems  to  perform  FI, 
there  are  two  alternative  subsystems,  but  due  to  some 
limitation  peculiar  to  one  hull  type,  only  the  first 
subsystem  can  be  used  on  it.  For  the  second  hull  type, 
however,  there  are  now  two  feasible  alternatives,  one 
with  the  first  subsystem  and  one  with  the  second. 


Suppose  noM  the  second  subsystem  has  features  which 
cause  the  designer  to  want  to  explore  an  alternative 
with  two  of  them  installed.  Either  one  can  effectively 
perform  the  function  FI,  so  the  subsystems  are  in 
parallel.  This  relationship  is  shown  in  Figure  (2-3). 


FIGURE  (2-3) 

SERIES  PARALLEL  MISSION  STRUCTURE  INCLUDING  SUBSYSTEMS 

There  are  now  four  alternative  designs  based  around 
function  FI,  the  hull  types  and  the  subsystems.  For 
the  first  hull  type  there  is  one  possible  subsystem, 
and  for  the  second  hull  type  there  are  three  possible 
combinations.  As  subsystems  are  selected  for  the  other 
functions  the  list  of  alternatives  will  grow.  Upon 
completion  of  this  process  for  every  mission  and 
function,  the  designer  has  a  set  of  alternative  ship 
system  specifications  or  lists  of  subsystems. 

SYNTHESIS  OF  SHIP  DESIGN  ALTERNATIVES 

The  lists  of  subsystems  are  not  yet  descriptions  of 
alternative  ship  designs.  A  great  deal  of  information 
must  yet  be  calculated  from  the  lists.  All  of  the 
subsystems  have  weight  and  take  up  volume.  In  addition 
many  have  electrical  power,  heating,  cooling  and 


ventilating  requirements.  The  size  of  the  support 
systems  must  be  calculated  from  these  requirements,  and 
the  size  of  the  ship  to  hold  all  this  plus  the  propul¬ 
sion  plant,  fuel,  food  and  stores  must  be  estimated 
before  an  alternative  ship  design  is  defined. 

Synthesis  models  are  series  of  calculations  which 
calculate  the  characteristics  and  dimensions  of  a  ship 
from  the  performance  requirements  such  as  maximum  top 
speed,  endurance  range  and  speed,  type  of  propulsion 
and  electrical  plants,  and  the  list  of  payload  items. 

A  synthesis  model  is  the  first  iteration  on  the  design. 
A  ship  synthesis  model  should  produce  enough 
information  about  the  ship's  physical  characteristics 
to  allow  the  designer  to  determine,  either  through 
direct  comparison  of  ship  characteristics  or  through 
further  analysis  of  those  characteristics  whether  that 
design  is  feasible  and  whether  it  is  better  or  worse  in 
comparison  to  other  alternatives. 

For  the  estimation  of  relative  mission 
effectiveness  the  synthesis  model  must  produce  enough 
information  about  the  ship's  hull  to  develop  a  hull 
form.  At  a  minimum  this  will  include  displacement  and 
overall  dimensions  of  the  hull.  Depending  on  the  hull 
type  and  method  of  developing  the  hull  shape,  other 
dimensions  or  coefficients  of  form  may  be  required  and 
must  be  produced  by  the  synthesis  model. 
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DEVELOPMENT  OF  A  HULL  SHAPE 


The  response  o-f  a  ship  to  the  ocean  waves  depends 
upon  the  shape  of  the  hul 1 . 

The  final  hull  form  will  be  a  compromise 
between  optimizing  for  arrangements,  minimum  resistance 
at  top  and  cruise  speeds,  seakeeping,  stability,  and 
other  limitations  which  may  arise  such  as  a  maximum 
limit  on  draft  for  various  ports.  None  of  this  can 
be  accurately  predicted  at  this  early  stage,  so  a 
"prototype"  hull  form  should  be  used  for  the 
calculation  of  ship  motions.  This  hull  will  not  have 
the  same  shape  as  the  final  hull,  nor  will  its  response 
be  identical,  either.  The  primary  criterion  is  that 
the  hull  shape  conform  to  physical  hull  parameters 
produced  by  the  synthesis  model,  and  that  there  be 
enough  information  about  the  shape  to  permit  the 
calculation  of  the  ship's  motions. 

ASSIGNMENT  OF  SUBSYSTEMS  TO  LOCATIONS 

The  amplitude  of  a  ship's  motions  varies  with 
location  on  the  ship,  the  largest  amplitudes  being  at 
the  bow  and  stern  and  the  smallest  near  the  center  of 
gravity  of  the  ship.  In  order  to  predict  the  motions 
experienced  by  a  subsystem,  its  location  on  the  ship 
must  be  known.  There  are  several  problems,  however, 
with  fixing  exactly  the  locations  of  subsystems. 
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First,  the  locations  will  almost  certainly  change  due 
to  revisions  and  compromises  necessary  for  a  variety  of 
reasons  imhich  will  only  come  to  light  later  in  the 
design  process.  Second  computing  the  motion  for  every 
location  would  involve  a  very  large  computational 
effort.  Third,  the  motions  at  a  location  will  not  be 
identical  to  the  final  design  motions  because  the  hull 
shape  being  used  is  not  the  same  as  the  final  hull 
shape,  as  stated  previously.  There  must  be  a  balance 
between  making  the  computations  reasonable,  and 
accuracy.  Only  a  few  representative  locations  should 
be  chosen.  Equipment  and  subsystems  are  assigned  to 
those  locations  that  they  will  be  closest  to,  and  once 
the  motions  are  computed,  the  motion  for  the 
appropriate  station  is  used  to  compute  performance 
degredation  for  a  subsystem. 

Some  subsystems  are  constrained  to  particular 
locations.  Hull  mounted  sonars  are  placed  at  or  near 
the  bow  to  be  as  far  away  from  propeller  and  machinery 
noise  as  possible.  Towed  sonars  must  be  deployed  from 
the  stern  to  avoid  entanglement  in  the  propeller. 
Helicopters  must  land  aft  of  the  deckhouse  and  masts. 
Radar  and  communication  antennas  are  generally  placed 
amidships  for  complete  coverage  of  the  horizon.  Weapons 
launchers  go  wherever  they  can  give  the  widest 
coverage.  These  constraints  must  be  kept  in  mind  when 


assigning  locations. 

CALCULATION  OF  SHIP  MOTIONS 

Current  techniques  are  all  numerical,  and  are  capable 
o-f  computing  a  ship’s  motions  for  only  one  combination 
of  speed,  sea  state  and  relative  wave  heading  at  one 
time.  Consequently,  motion  calculations  are  repetitive 
and  extremely  time  consuming  to  perform  by  hand.  Any 
program  which  allows  specification  of  speed,  sea  state, 
hull  form  and  location  for  computation  of  motion  and 
computes  the  necessary  motions  information  required  by 
the  subsystems’  performance  degredation  functions  may 
be  used  for  ship  motion  calculations. 

CALCULATION  OF  RELATIVE  MISSION  EFFECTIVENESS. 

At  this  point  the  designer  has  the  series-parallel 
structure  of  the  mission  relating  the  subsystems  to 
functions  and  then  to  the  mission,  and  motions  of 
locations  near  where  the  subsystems  are  expected  to  be 
for  each  speed,  sea  state,  and  wave  heading 
combination,  for  each  of  these  combinations,  one  after 
the  other,  the  relative  performance  of  each  subsystem 
is  determined  from  its  performance  degredation 
function.  The  mission  effectiveness  is  calculated  using 
equations  (2-1)  and  (2—2). 

The  designer  now  has  a  set  of  relative  mission 
effectivenesses  for  a  set  of  sea  conditions  for  a 
design  alternative.  The  process  is  repeated  for  all 
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the  alternatives,  and  the  seakeeping  results  become 
part  o-f  the  information  to  be  assessed  when  deciding 
which  alternative  is  the  "best.” 
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CHAPTER  4 


ILLUSTRATION  OF  THE  PROCESS  FOR  A  SINGLE  MONOHULL 

DESIGN 

To  illustrate  this  approach,  several  existing 
programs  were  used  and  one  program  was  written.  The 
existing  programs  are  as  follows:  A  ship  synthesis 
program  to  compute  ship  characteristics  based  on 
performance  requirements  and  a  given  set  of  payloads. 
Two  programs  are  used  which  calculate  a  “typical"  hull 
■form,  and  two  programs  which  use  the  hull  form  program 
output  and  calculate  the  ship's  response  to  long 
crested  head  seas.  The  program  which  was  written 
performed  the  calculations  for  mission  effectiveness. 
The  output  is  a  value  for  relative  effectiveness  of  a 
mission  area  at  a  given  speed  in  a  given  sea  state  for 
head  seas . 

MISSION  STRUCTURE 

This  illustration,  which  is  described  further  in 
Appendix  1,  used  a  simplified  set  of  mission 
requirements.  The  mission  requirements  document  which 
es  reproduced  in  Appendix  I  specified  warfare  and 
supporting  mission  areas,  and  physical  and  performance 


constraints  on  the  design.  Each  mission  area  Mas 
broken  down  into  major  functions  and  its  series  — 
parallel  mission  structure  developed  by  folloMing  the 
guidelines  in  chapter  2.  There  is  no  automated 
procedure  for  developing  this  mission  structure.  The 
mission  structure  was  developed  manually  by  the  author. 

Subsystems  Mere  next  selected.  Three  criteria  Mere 
used  by  the  author  in  choosing  the  subsystems.  First, 
that  they  performed  the  function,  second  that  they  Mere 
likely  to  meet  the  other  non— seakeeping  requirements, 
and  third,  they  Mere  in  the  data  base  for  the  synthesis 
model  that  Mas  used.  The  subsystems  Mere  then 
substituted  into  the  functions  and  the  mission 
structure  expanded,  and,  in  some  cases  revised.  It  Mas 
found  that  certain  subsystems  performed  more  than  one 
function.  When  this  happened,  the  tMO  separate 
functions  Mere  combined,  and  performed  by  one 
subsystem.  The  mission  diagrams  can  be  found  in 
Appendix  I,  along  Mith  the  complete  equipment  list. 
SYNTHESIS  OF  THE  DESIGN 

For  conventional  monohull  ships,  there  are  a  number 
of  computer  programs  Mhich  perform  the  synthesis 
calculations.  The  one  readily  available  to  the  author 
Mas  the  REED  model  ClOl  and  it  Mas  selected  for  use  in 
this  illustration  of  the  general  procedure.  The  REED 
model  is  a  regression  based  program  Mhich  meets  the 
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necessary  requirements  o-f  a  ship  synthesis  model;  it 
produces  the  ship’s  hull  -form  characteristics  necessary 
-For  the  development  o-F  a  hull  form,  plus  information 
which  allowed  the  author  to  compare  that  alternative  to 
the  other  non-seakeeping  constraints  specified  by  the 
mission  requirements. 

SHIP  MOTIONS 

For  a  conventional  monohull,  the  method  selected 
for  producing  a  prototype  is  the  use  of  the  two 
programs,  PREHULL  and  HUL6EN.  The  inputs  are  length, 
beam,  draft,  prismatic  and  sectional  area  coefficients 
and  depth  of  hull  at  midships.  The  result  is  a  hull 
form  which  is  not  optimized  in  any  particular  way,  but 
which  conforms  to  standard  US  Navy  design  practices. 
This  shape  would  be  the  starting  point  from  which  a 
designer  would  develop  a  hull. 

HULGEN  is  capable  of  producing  different  files  with 
different  types  of  information  about  a  hull.  One  of 
those  files  contains  information  on  the  sectional 
areas,  design  waterline  and  hydrostatic  properties. 

This  is  the  file  used  by  POSTHULL  to  create  the 
necessary  input  file  for  M0TI0N2D,  both  of  which  are 
further  described  below. 

The  locations  chosen  for  the  calculation  of  motions 
were  station  three  for  slamming,  stations  five, eight, 
fifteen  and  twenty  for  motions  of  displacement. 


velocity  and  acceleration.  Station  five  was  chosen  for 
ecfuipeent  on  the  forecastle.  Station  eight  was  chosen 
because  many  sensors,  the  bridge  and  CIC  are  often 
forward  of  amidships.  Slightly  aft  of  midships  are 
usually  located  main  propulsion  machinery  which  is 
largely  insensitive  to  motion.  Station  twenty  was 
chosen  for  the  deployment  of  towed  sonars  and  for 
landing  helicopters,  and  station  fifteen  was  chosen  for 
equipment  which  may  be  in  that  vicinity.  These 
locations  are  not  necessarily  the  most  representative 
choices,  nor  is  five  locations  necessarily  the  best 
compromise  between  computation  time  and  accuracy. 

POSTHULL  uses  the  output  file  from  HUL6EN  as  input 
and  produces  a  properly  formatted  input  file  for 
nOTION20  Cl  13.  The  program  uses  a  Lewis— form 
representation  C123  which  requires  draft,  maximum  width 
and  sectional  area  of  each  station.  It  computes  the 
transfer  function  for  heave  and  pitch  by  calculating 
the  response  of  the  ship  to  a  series  of  waves  of 
specified  frequency  and  unit  amplitude.  It  then 
computes  the  ship  response  spectrum  for  a  given  input 
spectrum,  the  vertical  motions,  velocity  and 
acceleration  spectra  and  statistics  at  specified 
locations  along  the  length  of  the  ship.  The  outputs  are 
the  response  spectrums  and  the  statistics  of  RHS,  one 
third  highest,  and  one  tenth  highest  amplitude  of 
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displacement,  velocity  and  acceleration  -For  each 
location,  and  expected  frequency  of  slamming  per  hour. 
MISSION  EFFECTIVENESS  CALCULATION 

A  program  was  written  which  implemented  the 
procedure  described  in  chapter  2.  The  program  was 
written  specifically  to  demonstrate  this  process  and 
consequently  does  not  have  all  the  features  the  author 
considers  necessary  for  regular  use.  It  is  included 
for  completeness  in  Appendix  III. 

Each  subsystem’s  location  was  input  along  with  the 
type  of  motion  which  caused  its  performance  to 
deteriorate.  The  appropriate  motion  for  that  location 
was  found  and  used  to  compute  the  performance 
degredation  from  the  subsystem’s  performance 
degredation  function.  The  process  for  computing  the 
mission  effectiveness  described  in  chapter  2  was 
implemented  and  the  mission  effectiveness  for  one  sea 
state  and  two  speeds  calculated.  The  numerical  results 
are  in  appendix  I. 
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CHAPTER  5 


CONCLUSIONS  AND  RECOMMENDATIONS 


CONCLUSIONS 

This  illustration  o-f  the  approach  has  demonstrated 
that  it  is  possible  to  carry  out  the  modified  design 
process  proposed  in  chapter  3.  Calculations  can  be 
performed  in  a  reasonable  time  through  the  use  of 
computer  programs  to  carry  out  the  numerical 
computations. 

The  actual  mission  effectiveness  procedure  is  not 
directly  linked  to  any  synthesis,  hullform  generation, 
or  ship’s  motions  calculation  procedures.  The  purpose 
of  those  three  steps  is  to  compute  a  set  of  motions  for 
a  ship  which  will  carry  the  subsystems.  However,  for 
the  complete  evaluation  of  an  alternative  on  both 
seakeeping  and  non  seakeeping  criteria  a  synthesis 
model  of  some  form  is  essential. 

The  value  produced  by  this  approach  is  more  of  a 
figure  of  merit  than  a  true  measure  or  the  ability  of  a 
ship  to  perform  a  mission  in  a  given  set  of  sea 
conditions  relative  to  calm  water.  Because  the  hull 
and  the  locations  of  subsystems  cannot  be  expected  to 


be  identical  to  the  -final  design,  the  actual  mission 
e^-fecti  veness  cannot  be  expected  to  be  the  same.  For 
two  alternative  designs  which  compare  closely,  a  small 
di-fference  in  value  would  be  meaningless. 

The  lack  o-f  -firm  performance  degredation  infor — 
(nation  on  subsystems  severely  limits  the  usefulness  of 
this  approach.  Until  such  time  at  there  is  more 
accurate  data  on  subsystems,  it  is  the  author's  opinion 
that  the  approach  be  used  with  great  caution,  and 
should  not  be  relied  upon  solely  when  evaluating 
alternative  designs. 

RECOMMENOAT I ONS 

Efforts  are  underway  to  collect  better  information 
on  the  performance  degredation  of  subsystems.  This 
work  should  be  continued. 

The  mission  structure  is,  to  a  great  extent, 
independent  of  scenario.  In  different  scenarios  for  a 
mission,  different  subsystems  may  be  used,  and  one  or 
more  may  not  be  used  at  all.  While  placing  subsystems 
in  parallel  implies  that  the  performance  of  all  of  them 
is  not  necessary  to  complete  a  mission,  the  underlying 
assumption  is  that  any  of  the  parallel  subsystems  is 
equally  capable  of  performing  its  function.  This  may 
not  be  true.  The  effect  of  this  assumption  on  the 
mission  effectiveness  needs  to  be  studied  carefully. 

Further  work  on  implementing  this  procedure  with 


other  types  of  hull  shapes  should  be  carried  out  to 
provide  the  designer  with  the  tools  to  compare 
different  alternatives. 

This  procedure  should  be  implemented  for  use  with 
the  design  tools  already  in  place  at  NAVSEA. 
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APPENDIX  I 


EXAMPLE  OF  THE  APPROACH 


MISSION  REQUIREMENTS 

In  order  to  demonstrate  hoM  the  mission 
eT-f ecti veness  calculation  process  can  be  applied  to 
feasibility  studies,  a  highly  simplified  feasibility 
study  is  described  here.  It  begins  with  a  set  of 
mission  requirements  and  constraints  which  are  taken 
from  an  example  given  in  a  Combat  Systems  Design 
course,  taught  in  the  summer  of  1983.  The  directions 
are  reproduced  here. 

DESIGN  DIRECTION  FOR  “LOW  MIX”  SURFACE  COMBATANT 
1.  Mission  Statement 

The  "low  mix"  Surface  Combatant  is  to  be  a  replacement 
for  the  FF  1040,  FFG  1,  FF  1052  Class  frigates.  This 
ship’s  primary  mission  will  be  to  escort  low  value 
naval  task  forces  (replenishment  and  amphibious)  and 
mercantile  fleets  in  relatively  low— threat 
environments.  The  primary  threat  is  from  diesel  and 
nuclear  submarines. 
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2.  Per-formance  Requirements 

a.  ASW  -  Defend  task  force  against  enemy  attack 
submarines. 

b.  AAW  —  Limited  defense  for  task  force.  Point 
defense  for  own  ship. 

c.  SUM  -  Defend  task  force  against  surface  forces. 
Shore  bombardment. 

d.  C^I  —  assist  in  control  of  task  force/mercantile 
force  assuming  presence  of  at  least  one  mid  mix 
ship  serving  as  force  commander. 

e.  Mobility  —  Maximum  sustained  speed  —  28  knots 
(baseline).  Minimum  range  4,500  NM  at  18  knots. 

f.  UNREP  -  Refuel  and  receive  stores  and  munitions 
while  underway. 

3.  Constraints 

a.  Acquisition  Cost  —  This  ship  will  be  needed  in 

large  numbers  to  replace 
retiring  assets.  A  follow 
ship  acquisition  cost 
constraint  of  ♦300M  (FY  80) 
is  established. 

b.  Ship  Displacement  —  full  load  displacement  will 

not  exceed  4000  tons. 

c.  Manning  -  Number  of  accommodations  will  not  exceed 

225. 
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Performance  requirements  a  through  f  above  are 
mission  areas.  Three  mission  areas  are  considered 


supporting  mission  areas;  they  are  Mobility  (MOB), 
Underway  Replenishment  (UNREP)  and  Command,  Control, 
Communications,  and  Intelligence  gathering  (C^I).  All 
three  warfare  areas  are  diagrammed  in  parallel.  As 
long  as  the  ship  is  capable  of  performing  one  mission, 
it  is  still  capable  of  performing  some  part  of  its 
overall  mission  as  escort.  This  is  an  arbitrary 
decision  made  primarily  to  simplify  this  illustration 
of  the  process.  It  is  quite  possible  that  some  other 
emphasis  (ie  ASW  is  essential  to  perform  the  entire 
mission)  on  warfare  areas  would  be  equally  correct. 


replenish  while  underway,  or  with  no  communications 
capability.  The  ship  can  still  perform  its  mission, 
although  in  a  degraded  capacity  if  any  one  or  even  two 
of  the  warfare  mission  areas  could  not  be  performed.  A 
rather  critical  assumption  is  that  a  particular  warfare 
area,  say  Anti-Submarine  Warfare  (ASW>  will  not  always 
be  necessary.  The  choice  of  which  functions  must  be 
performed  is  a  matter  of  Judgement  on  the  part  of  the 
designer.  Two  equally  acceptable  alternatives  to  the 
parallel  relation  are:  1)  that  all  three  warfare 
mission  areas  would  be  in  series.  2)  since  the  primary 
threat  is  submarines,  ASW  is  in  series  and  that 
Anti— Air  Warfare  (AAW)  and  Surface  Warfare  (SUW)  are  in 
paral lei . 

There  are  five  basic  functions  common  to  all  three 
warfare  areas  (figure  A-2)  detection,  identification, 
tracking,  engaging  or  attacking,  and  evaluating  the 
attack  for  success.  In  each  warfare  area  these 
functions  are  performed  different  subsystems.  All  five 
functions  are  considered  vital;  no  one  function  can 
fail  completely  without  the  warfare  capability 


f 

1/C.  1  1 

1 

DETECT^ - ^ IDENTIFY^ - ^ TRACKS -  EN6ABEI - 1  EVALUATE 


FIGURE  (A-2) 

DIAGRAM  OF  FUNCTIONS  NECESSARY  TO  EACH  WARFARE  AREA 
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failing  completely. 

SELECTION  OF  SUBSYSTEMS  FROM  TECHNOLOGY  BASE 

At  this  point  it  becomes  necessary  to  use  the 
available  technology  base,  the  data  base  for  the  Reed 
ship  synthesis  model.  In  conjunction  with  developing 
the  functional  relationships,  selection  of  specific 
subsystems  for  the  synthesis  model  are  done  for  one 
alternative.  The  complete  list  of  subsystems 
for  the  warfare  mission  areas  is  shown  in  table  (A-1). 
The  list  of  subsystems,  along  with  the  performance 
parameters  were  input  to  the  REED  model.  The  resulting 
REED  model  output  is  shown  in  appendix  II. 

The  subsystems  were  incorporated  into  separate 
series-parallel  functional  diagrams  for  each  mission. 
The  resulting  diagrams  are  shown  in  figures  (A— 3) 
through  (A-5> .  The  relationships  shown  are  not 
necessarily  the  only  passible  relationships. 

Some  of  the  subsystems  depicted  in  these  diagrams 
are  actually  functions.  An  example  is  launching  and 
recovering  a  helicopter.  This  involves  the  hull,  on 
which  the  helicopter  lands,  support  personnel  on  the 
ship,  the  helicopter,  and  its  pilot.  The  operation  is 
generally  viewed  as  a  single  function  or  subsystem,  and 
performance  degredation  is  estimated  for  the  entire 
operation,  not  for  the  many  steps  involved  C51.  The 
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FUNCTION 


SUBSYSTEM 


ISW 

DETECT,  IDENTIFY 
TRACK,  EVALUATE 

ENGAGE 

^AW 

DETECT,  TRACK 

EVALUATE 

IDENTIFY 

ENGAGE 


>UU 

DETECT,  TRACK 
TRACK,  EVALUATE 

ENGAGE 


SQS-56  sonar 
SQR-19  sonar 
2  LAMPS  helos 
LAMF>S 
t1K-Z2  T.T. 

SPS-49  Radar 
SPS-10  Radar 
IFF  system 
SLQ-32 
NATO  SEA 
SPARROW 
MK-86  GFCS 
MK-45  naval  gun 
CIWS 

SPS— 10  radar 
LAMPS  helo 
MK-86  GFCS 
MK— 45  naval  gun 
HARPOON  ASCM 
TABLE  (A-1) 


SUBSYSTEMS  FOR  THE  DESIGN  ALTERNATIVE 


"subsystefli"  designated  PERS,  for  personnel,  is  included 
where  the  subsystem  is  operated  by  personnel,  and  their 
performance  can  affect  the  function  that  the  subsystem 
and  personnel  together  perform.  The  signals  processed 
by  the  sonar  are  evaluated  by  the  sonar  operators.  That 
information  is  passed  on  to  other  personnel  who 
combine  it  with  information  from  other  sources  in 
identifying  and  tracking  a  target,  and  controlling  an 
attack. 


ASM 


(H 

SQS-561 — 1 

PERS 

|-|SQR-19| — 

M 

1  LAUNCH/LAND  1 — 

jPERSl 

H 

DETECT,  TRACK,  IDEhlTIFY, 


_ /  V. 

EVALUATE 


ENGAGE 


j 


FIGURE  (A-3) 

ASW  MISSION  EXPANDED  INTO  ITS  COMPONENT  SUBSYSTEMS 
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AAW 


r— 

MK-Q6 

MK-45 

— J  PERS  - 

GFCS 

LWG 

-  PERSi— 

-fcIwS 

4slq-32MpersI— 1 

— 

-  NATO 

SEA  SPARROW  ^ - 1 

-'v- - - __y 


DETECT,  TRACK,  IDENTIFY  ENGAGE 

EVALUATE 


FIGURE  (A-4) 

AAM  MISSION  EXPANDED  INTO  ITS  COMPONENT  SUBSYSTEMS 


SUM 


DETECT,  TRACK,  IDENTIFY,  EVALUATE  ENGAGE 


FIGURE  (A-5) 

SUM  MISSION  EXPANDED  INTO  ITS  COMPONENT  SUBSYSTEMS 
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Two  programs  are  used  to  generate  a  hull  -form.  The 
first  is  PREHULL  which  accepts  the  basic  hull 
parameters  of  length,  beam  draft,  depth  of  hull  at 
station  ten,  and  the  prismatic  and  midship  section 
coeffients,  and  creates  and  input  file  for  HULGEN, 
which  is  an  interactive  hull  form  generation  program. 
Both  of  these  (as  well  as  POSTHULL  and  nOTION2D,  used 
in  the  motion  calculations)  are  part  of  the  NAVSEA  Hull 
Form  Design  System  currently  used  for  preliminary  hull 
design.  The  file  produced  by  PREHULL  describes  a 
conventional  monohull  which  conforms  to  standard  Navy 
design  practices.  The  output  of  the  HULGEN  programs 
forms  the  input  to  the  programs  which  compute  ship's 
motions. 

five  locations  for  calculating  motions  for 
subsystems  were  chosen  for  this  demonstration.  The 
locations  were:  Station  five,  one  quarter  of  the 
ship’s  length  from  the  bow  for  equipment  on  the 
forecastle.  Station  eight  for  combat  systems,  CIC 
equipment,  and  equipment  on  masts,  etc.  Stations 
fifteen  and  twenty  for  equipment  aft,  particularly 
helicopter  launching  and  landing,  and  towed  sonar 
systems. 
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SYSTEM/FUNCTION 


STATION 


SQS-56 

(none  subject 

SQR-19 

20 

LAMPS  hel i copter 

15 

MK— 32  Torpedo  Tubes 

8 

SPS-49 

8 

SPS-IO 

8 

SLQ-32 

8 

MK-8A  GFCS 

8 

NATO  SEA  SPARROW 

20 

MK-45  LWG 

5 

CIWS 

15 

HARPOON  ASCM 

5 

UNREP 

8 

TABLE  (A-2) 

LOCATION  OF  MAJOR  SUBSYSTEMS 


Far  this  example  the  MOTION  2D  Cl 13  program  was 
used  to  compute  ship  response  to  waves.  The  program  is 
limited  to  long  crested  head  seas,  and  the  Pierson 
Moskowitz  sea  spectrum  -formula  was  used  to  generate 
spectra. 

CALCULATION  OF  MISSION  EFFECTIVENESS 

Up  to  this  point,  all  three  missions  have  been 
dealt  with  in  order  to  -fully  illustrate  the  process. 

The  process  o-f  calculating  mission  effectiveness  is 


identical  -For  each  mission  so  only  one,  the  ASW  mission 
calculations  will  be  shown  here. 

Normally,  the  designer  would  examine  all  sea  state 
and  speed  conditions  considered  to  be  of  interest  or 
importance  for  comparison.  Only  two  sample 
calculations  were  performed  for  sea  state  six,  fully 
developed,  long  crested  head  seas,  at  five  and  fifteen 
knots.  The  procedure  for  both  is  identical,  as  is  the 
procedure  for  any  other  set  of  sea  states  and  speeds. 

A  program  was  written  to  demonstrate  the 
calculation  process.  The  program  has  three  inputs  to 
it.  First,  a  list  of  the  systems,  whether  they  are  in 
series  or  in  parallel,  and  their  assigned  locations. 
Second,  the  vertical  motions,  velocities  and 
accelerations  of  the  ship  at  each  of  four  stations  plus 
expected  frequency  of  slam.  Third,  performance 
degredation  functions  for  the  functions  and  subsystems. 

It  must  be  recognized  that  any  subsystem  not 
affected  by  motion  need  not  be  included  in  the  diagram 
for  the  purpose  of  calculating  mission  effectiveness. 
The  effects  of  motion  on  command  control  equipment, 
mainly  communications  gear,  navigational  electronics, 
and  computers  are  negligible.  The  effects  on  the 
personnel  who  operate  that  equipment  and  evaluate 
information  coming  from  it  cannot  be  ignored,  however. 
Results  of  the  calculation  are  shown  in  figures  (A-6) 
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and  (A--7) .  They  are  summarized  below 


Mi ssion 
Sea  state 
Speed  (kts) 
Effectiveness 


ASW 


h  6 

5  15 

1.0  0.58 


The  results  are  interpreted  as  follows.  In  sea  state 
six,  for  a  speed  of  five  knots,  the  ship  may  perform 
the  ASW  mission  about  as  well  as  in  calm  water.  In  sea 
state  six,  for  a  speed  of  fifteen  knots,  the  ship  may 
see  a  reduction  in  effectiveness  of  about  half.  The 
computer  program  is  capable  of  producing  seven  or  more 
digits  of  accuracy.  On  the  other  hand,  the  data  and 
the  assumptions  involved  with  the  procedure  reduce  the 
accuracy  greatly. 


FlC  (A-6) 

MISSION  EFFECTlYENesS,  FIVE  f^TS.  SEA  STATE  S/X 
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THE  EFFECTIVENESS  RELATIVE  TO  THIS  FEASIBLE 
DESIGN  CALM  WATER  PERFORMANCE  IS 


FIG  (A- 7) 

mis5;on  effectiveness,  fifteen  kts.  sea  state  six 
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APPENDIX  II 


REED  MODEL  OUTPUT 
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REED  MODEL  OUTPUT  (CONTD) 
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MISSION  EFFECTIVENESS  CALCULATION  PROGRAM 
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PR06RAN  PERFORH 
C 

CHflRACTER»64  SOURCE,  ROTNS,  DBLIST,  DBASE,  REPLY 
CHARACTER»12  NAI1ES(50,2) .DBLlSdOO) 

CHARACTER#!  TYPE(50,2) 

1HTE5ER«2  RONS,  ISCURC,  IHOTNS,  IDBLIS,  {DBASE 
1 NTE6ER*2  LOCATE (50,2!, DBADOR ( 1 001 , DBS! ZE 
REAL  VALUES (50) 

C 

COnnON  /INOUT/  IS0URC,I)10TNS,IDBLIS,IDBAS£ 

CORHON  /3YSTEI1/  NAMES,  TYPE, LOCATE,  VALUES 
COnnON  /SYSBAT/DBSIZE,DBL{S,DfiADDR 
C 

C  ENSURE  THAT  CHARACTER  VARIABLES  ARE  CLEARED  AND  HAVE  ONLY  BLANKS 
C 

CALL  EMPTY  (SOURCE) 

CALL  EMPTY  (MOTNS) 

CALL  EMPTY  (DBLIST) 

CALL  EMPTY  (DBASE) 

CALL  EMPTY  (REPLY) 

C 

HRITE  (#.50) 

50  FORMAT  ill,'  BEFORE  ENTERINS  THIS  PR06RAM  YOU  SHOULD  HAVE  HANDY’, 
4  ’  THE  FILE’,/,’  NAMES  FOR  THE  DATA  BASE,  SHIP  MOTIONS  AND  THE’, 

4  /,’  SYSTEMS  DIAGRAM,  OR  THE  DIAGRAM  IN  HAND.’) 

C 

HRITE  (#,55) 

55  FORMAT  (//,’  NILL  YOU  ENTER  SYSTEM  INFORMATION  FROM  TERMINAL  ’, 

4  ’OR  FILE’,/,’  TYPE  IN  NHICH’) 

READ  i»,60)  REPLY 
60  FORMAT  (A64) 

C 

IF  (REPLY  ,EQ.  ’FILE’)  THEN 
4RITE  (»,65) 

65  FORMAT  (/,’  NHAT  I:  THE  FILE  NAME?  ’! 

READ  :t,60)  SOURCE 
ISOURC  =  10 

OPEN  !  ISOURC. FILE=SOURCE,STATUS=’CLD’) 

READ  (ISCURC, 70) RONS 
70  FORMAT  1BN,I10) 

ELSE 

NfilTE  (#,75) 

75  FORMAT  (/,’  ENTER  THE  TOTAL  NUMBER  OF  ELEMENTS  IN  THE  SYSTEM’! 
READ  (#,70)  RONS 
END  IF 
C 

NRITE  (#,e0) 

30  FORMAT  (//,’  NHAT  FILE  IS  THE  SHIP  MOTIONS  STORED  IN’’  ’) 

READ  (#,60)  MOTNS 
i MOTNS  =  20 

OPEN  ! IMOTNS.FILE=HOTNS,STATUS=’OLD’ ) 

r 


bl 


C  tttltttltttttltltttltttltttttilttttttlltitMttUIttllllllttlllttttll 

C 

SUBROUTINE  MINSUBIRONS, REPLY) 

C 

CHARACTER»64  REPLY 

CHARACTERt 1 2  NAHES ( 50 , 2) , D8L IS ( 1 00 ) 

CHARACTER#!  TYPE (50, 2) 

INTE6ER»2  L0CATE(50, 2), RONS, ISOURC.INOTNS.IDBLIS, (DBASE 
INTE8ER*2  DBADDR(100),SBSIZE 
REAL  VALUES (50),  EFFECTIVENESS 
C 

COnnON  /INOUT/  ISCURC,inOTNS,IDBLIS,ID6ASE 
COmON  /SYSTEM/  NAHES,  TYPE, LOCATE, VALUES 
COnnON  /SYSDAT/DBSIZE.SBLIS,DBADDR 
C 

C  ENSURE  THAT  CRITICAL  ARRAYS  START  OUT  ZEROED  OUT 
C 

CALL  ZERO  (RONS) 

C 

C  &ET  THE  SYSTEM  LIST  AND  RELATIONSHIPS  EITHER  FROM  FILE  OR  FROM 
C  THE  TERMINAL.  STORE  THE  FILE  IF  DESIRED 
C 

CALL  FILNAMS  (ROHS,REPLY) 

C 

C  READ  IN  THE  MOTIONS  OF  THE  SHIP,  CALCULATE  THE  PERFORMANCE 
C  EFFECTIVENESS  OF  THE  SEPARATE  EQUIPMENTS  OR  SUBSYSTEMS 
C 

CALL  SETHOT  (RONS) 

C 

C  PERFORM  THE  CALCULATION  FOR  MISSION  EFFECTIVENESS 
C 

CALL  MAINCALC  (RONS,  EFFECTIVENESS) 

r 

w 

C  NRITE  THE  OUTPUT  TO  A  FILE 
C 

CALL  OUTPUT  (RONS,  EFFECTIVENESS) 

C 

RETURN 

END 

C 
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C  tttltttltttltlMttttttltlidtlHtllitlttllHttHlllltltlllttltllllt 

C 

SUBROUTINE  EMPTY (STR) 

C 

CHflRACTER»64  STR 
C 

C  SET  STR  EQUAL  TO  BLANKS-THERE  ARE  64  BLANKS 
C 

STR=’ 

tr 

RETURN 

END 

C 

C 

SUBROUTINE  ZERO  IRONS) 

C 

CHARACTERfl2  NAnES(S0,2),EHPTY 
CHARACTER*!  TYPE (50, 2) 

INTEGER*2  LOCATE(SO,2),  RONS 
REAL  VALUESISOl 

r 

CCMGN  /SYSTEM/  NAMES, TYPE, LOCATE, VALUES 

r 

w 

DO  1010  1*1,50 
DO  1020  J*l,2 
NANESII.J)  »  ’ 

LOCATEII.J)  =  0 
1020  CONTINUE 

VALUES (1)  »  0.0 
1010  CONTINUE 
C 

RETURN 

END 

r 

u 
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C 

SUBROUTINE  FILNflHSIRONS, REPLY) 

C 

CHARACTERt64  REPLY, EMPTY 
CHARACTER»12  NAMES (50, 2), SOURCE 
CHARACTERil  TYPE(50,2) 

INTEGERt2  LOCATE(50,2>,RONS,I,ISOURCE,IMOTNS,1BBASE,IDBLIS 
C 

COMMON  /INOUT/  IS0URCE,IM0TNS,IDBLIS,1DBASE 
COMMON  /SYSTEM/  NAMES, TYPE.I.QCATE, VALUES 
C 

NRITE  (EMPTY, 101) 

101  F0RMAT(44)() 

SOURCE  =  EMPTY 
C 

IF  (REPLY  .EQ.  ’TERMINAL’)  THEN 
DO  2010  1=1, RONS 
10  NRITE  (♦,201) 

201  FORMAT  (//,’  BLOCK  NAME?  ’) 

READ  (♦.210)  NAMES(I,1) 

210  FORMAT  (A12) 

IF  til  .ST.  1).AMD.(NAHES(I-1,1>  .E8.  NAMESd,!)))  THEN 
TYPEd.l)  =  TYPEd-1,1) 

ELSE 

C 

NRITE  (♦,203) 

203  FORMAT  (/,’  IS  THIS  A  SERIES  OR  PARALLEL  BLOCK’,/, 

I  ’  TYPE  S  OR  P  ’) 

READ  (♦,212)  TYPEd,l) 

212  FORMAT  (Al) 
r 

END  IF 

r 

-iSITE  -♦,:02) 

202  FORMAT  (/,’  ELEMENT  NAME?  ’) 

READ  (♦,210)  NAMESd, 2) 

NRITE(^,204) 

204  FORMAT!/,’  IS  THIS  ELEMENT  ANOTHER  BLOCK,  OR  EBUIPMENT/’, 

(  ’SUBSYSTEM’,/,’  TYPE  B  (BLOCK)  OR  E  (EBUIPMENT)’) 

READ  (♦,212)  TYPE  d,2) 

NRITE  (♦,205) 

IF  (TYPE  (1,2)  .£3.  ’E’)  THEN 

205  FCRMAT(/,  ’  FOR  APPROXIMATE  LOCATION,  TYPE  THE  MENU  ’,/, 

4  ’  NUMBER  YOU  THINK  THE  EBUIPMENT/SYSTEM  IS  NEAREST’) 

NRITE  (♦,206) 

206  FORMAT!/, 5X,’  1’, 5X, ’NEAR  BON’,/, 5X,’  2’, 5X,’NEAR  MIDSHIPS’, 

4  /,5X,’  3’,5X.’NEAR  STATION  I5’,/,5X,’  4’,5X,’AT  STERN’) 

READ  (♦,213)  LOCATEd.l) 

213  FORMAT  (BN, 13) 

ELSE 

LOCATEd,!)  =  0 
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c 


END  IF 


20  NRITE  (»,207)  I 

207  FORNAT  !///,5)(, ’RON  ’,13,//,’  BLOCK’, 7!t,’NANE’,71(, 
i  ’SER/PAR  TYPE  LOCATION’) 

NRITE  (»,20B)  NAHES(I,l),NAHES«I,2),TyPE{I,l),TYPE(I,2l, 

&  LOCATEd,!) 

208  F0RNAT(1X,A12,A12,3X,A1,7](,A1,7X,I3} 

NRITE  (f,209) 

209  FORNATI/,’  IS  THE  ABOVE  CORRECT?  TYPE  Y  OR  N  ’) 

READ  (♦.214)  REPLY 

214  FORNAT(A) 

IF  (REPLY  .EQ.  ’N’)  THEN 
SO  TO  10 

ELSE  IF  (REPLY  .ED.  ’O’)  THEN 
60  TO  99 

ELSE  IF  (REPLY  .NE.  ’Y’)  THEN 
NRITE  (*,9011  REPLY 

901  FORNAT (/,’  YOUR  REPLY  ’,A12,  ’  NAS  NOT  UNDERSTOOD’) 

REPLY  -  EHPTY 
GO  TO  20 
END  IF 
2010  CONTINUE 

r 

3010  TYPE  ♦,  ’  ’ 

TYPE  »,  ’  DO  YOU  NANT  TO  SAVE  THIS  IN  A  FILE?’ 

TYPE  ♦,  ’  TYPE  Y  OR  N’ 

READ  (t,214)  ANS 
IF  (ANS  .£a.  ’Y’)  THEN 
TYPE  »,  ’  ’ 

TYPE  »,  ’  NHAT  IS  FILE  NAHE?’ 

READ  (♦,2141  SOURCE 
ISOURCE  =  10 

OPEN  (ISOURCE,  FILE  -  SOURCE,  STATUS  =  ’UNKNCNN’) 
j(RITE  (ISOURCE, 450)  RONS 
450  FORNAT (15) 

DO  1=1,  RONS 

NRITE  (ISOURCE, 350)  NAHES(I,1),NARES(I,2),TYPE(I,1), 
k  TYPEd, 2), LOCATEd,!) 

END  DO 

ELSE  IF  (ANS  .NE.  ’N’)  THEN 
TYPE  ♦,’  ’ 

TYPE  ♦,  ’  ANSNER  NOT  UNDERSTOOD’ 

60  TO  3010 
END  IF 
C 

ELSE  IF  (REPLY  .EQ.  ’FILE’)  THEN 
DO  2020  1=1, RONS 

READ  (IS0URCE.350)  NARESd,ll,NARESd,2),TYPE(l,l),TYPEd,2) 
i  .LOCATE!!, 1) 

350  F0RHAT(2A12.2(X,A1)I3) 

2020  CONTINUE 
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c 


END  IF 
BO  TO  999 


99  NRITE  (1,910) 

910  FORHAT  (////,’  ENDING  PROGRAH.’,//, 
STOP 

999  CONTINUE 
RETURN 
END 


QUIT’) 
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C 

SUBROUTINE  HAINCALC  (RONS,  EFFECTIVENESS) 

C 

CHARACTERtlZ  NAMES (50,2) 

CHARACTER*!  TYPE(50,2),0UHHY 
INTE6ER*2  PTR,  STPTR,  ISTACK(IO),  I,  LOCATE(5O,2),R0NS 
REAL  CURVAL,VALUES(50),RSTACK(10),  EFFECTIVENESS 
C 

CGHNON  /SYSTEM/  NAMES,  TYPE,  LOCATE,  VALUES 
C 

C  INITIALIZE  VARIABLES 

C 

C 

CURVAL  =  0.0 
PTR  =  1 
STPTR  =  1 
C 

C  START  CALCULATION  LOOP 
C 

C  FIND  END  OF  BLOCK 
C 

1  IF  ((PTR+1  .LE.  RONS)  .AND.  (NAMES(PTR*1,1)  .ED.  NAMES(PTR,1))) 
&  THEN 

PTR  »  PTR+1 
60  TO  1 
END  IF 
C 

C  IF  THIS  IS  A  BLOCK,  THEN  SET  CURRENT  VALUES  AND  STATUS  ON  THE 
C  STACK  AND  START  THE  PROCEDURE  ABAIN,  RECURSIVELY 
C 

.  2  IF  (VALUES (PTS)  .EQ.  0.0)  THEN 
IF  (TYPE(PTR,2)  .EQ.  ’B’)  THEN 
ISTACK (STPTR)  =  PTR 
RSThCK{STP’R)  =  CURVAL 
CURVAL  =0.0 
STPTR  =  STPTR  +1 

CALL  FINDBLOK  (R0NS,PTR,NAHES!PTR,2), NAMES) 

60  TO  I 
END  IF 
ELSE 

IF  !TYPE(PTR,1)  .EQ.  ’S’)  THEN 
IF  (NAf1ES(PTR+l,l)  .NE.  NAHES(PTR,1))  THEN 
CURVAL  =  VALUES (PTR) 

ELSE 

CURVAL  =  CURVAL  *  VALUES (PTR) 

END  IF 

ELSE  !F  (TYPE(PTR,1)  .EQ.  ’P’)  THEN 
IF  (NAMES (PTR+1,1)  .NE.  NANES(PTR,1))  THEN 
CURVAL  =  l.O-VALUES(PTR) 

ELSE 

CURVAL  =  CURVAL* (1.0-VALUES (PTR)) 


END  IF 
ELSE 

TYPE  ♦,  ’  ERROR  IN  SYSTEH  FUNCTION  LIST  AT  LINE  ’,PTR 
STOP 
END  IF 

IFKPTR  .6T.  1)  .AND.  (NAHESIPTR-l,!)  .EO.  NAnES(PTR,l)))THEN 
PTR  =  PTR-1 
SO  TO  2 

ELSE  IF  (PTR  .ST.  1)  THEN 
IF  (STPTR  .ST.  1)  THEN 

IF  (TYPE  (PTR.l)  .EQ.  ’PM  THEN 
CURVAL  =  1.0  -  CURVAL 
END  IF 

PTR  =  ISTACK(STPTR-l) 

VALUES (PTR)  =  CURVAL 
CURVAL  =  RSTACK(STPTR-l) 

STPTR  =  STPTR- 1 
SO  TO  2 
END  IF 
ELSE 

IF  (TYPE(PTR,1I  .EO.  ’P’)  THEN 
EFFECTIVENESS  =  I-CURVAL 
ELSE 

EFFECTIVENESS  =  CURVAL 
END  IF 
80  TO  999 
END  IF 
END  IF 
C 

999  CONTINUE 
RETURN 
END 


C 
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C  itHIttllttlltllttlttltttittttlttllMiltliHHitHlilttltttttltttt** 

C 

SUBROUTINE  FINOBLOK  (RONS, PTR.PTRVALUE, NAMES) 

C 

CHARACTERtl2  NAnES(50,2),  PTRVALUE 
1NTEBER*2  RONS,  PTR 
C 

PTR  =  1 

DO  NHILE  ((PTRVALUE  .NE.  NAMES (PTR, D). AND. (PTR  .LE.  RONS)) 

PTR  =  PTR  t  1 
END  DO 

IF  ((PTRVALUE  .NE.  NAHES(PTR,1)).AND.(PTR  .EQ.  RONS))  THEN 
STOP  ’ERROR  FROM  FINDBLQK  IN  NAMES’ 

END  IF 
RETURN 
END 


70 


c*>  «.**>  c*> 


c 

SUBROUTINE  SETNOT  (RONS) 

C 

CHARACTER*12  NAHES(50.2),DBLIS(100),  RECNANE 
CHARACTER*!  TYPE (SO, 2),  DURNY 
C 

INTE6ER«2  LOCATE (50, 2), RONS, ISQURC,  in0TNS,I08LIS.I,J,ir 
INTE6ER*2  L,DBADDR(100),RECNUN,  REC»OT,RECNT,flOTCOL 
!NTE6ER*2  (DBASE,  ERROR,  DBSIZE 

C 

REAL  VALUES(SO),  H0TNS(17),  RECDEG(10,2),  SYSHOT,  LINTERP 
REAL  HinOT,  LOnOT,  HIDES,  L0DE6,  nOTSPC(2S,13),  PRFSPC(2S) 

REAL  HSUfl,  PSUH,  ERRLIB,  VARIANCE 
C 
C 

COnnON  /INOUT/  isourc,irotns,idblis,idbase 
CONHCN  /SYSTEH/  NAMES. TYPE, LOCATE. VALUES 
COHHON  /SYSDAT/DBSIZE,DBLIS,DBAODR 
C 
C 

ASSI8N  999  TO  ERROR 
C 

C  READ  IN  NOTIONS  OF  LOCATIONS 
C 

READ  (INOTNS.lOO)  (NflTNS(J),JM,17) 

100  FORMAT (17F7. 3) 

C 

SO  IM,2S 

READ  (IMOTNS,  105)  (MOTSPCd.J),  3*1,13) 

105  F0RMAT(F10.4,12E14.7I 
END  SO 

READ  IN  DATABASE  LIST  AND  ADDRESSES 
SO  1=1, DBSIZE 

READ  (ISBLIS.llO)  DBLI3!I),SBADDR(I) 
no  F0RHAT(A12,I5) 

END  SO 
5  FORMAT (A) 

C 

C  CALCULATE  MOTION  DESREDATION 

r 

W 

SO  2030  1*1, RONS 
IF  (TYPE  (1,2)  .EQ.  ’EM  THEN 

C  FIND  ADDRESS  OF  ELEMENT 

3*1 

10  IF  (NAMESd.Z)  .EQ.  DBLIS(3))  THEM 
RECNUM  =  CBAD0R(3) 

60  TO  20 

ELSE  IF  (3  .EQ.  SSRONS)  THEN 
SO  TO  ERROR 
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ELSE 
J=J+1 
SO  TO  10 
END  IF 

C  RECORD  ADDRESS  FOUND  GET  NOTH  DEG. 

20  READ  (IDBASE,REC=RECNUH,FNT=1000>  RECNANE.RECnOT, (REC0EG(K,1), 
k  RECDE6(K,2},  K=l,10),  ERRLIfl 
1000  F0RnAT(A12,I2,10(F5.3,F7.2),F7.31 

C 

C  NON  HAVE  EQUIP  HOTIQHS  FUNCTION,  FILL  VALUES  NITH  I10T10HS 
C 

C  FIND  CORRECT  COLUHN  OF  NOTIONS 

C 

IF  (RECHOT  .ED.  1)  THEN 
HOTCOL  =  3 

ELSE  IF  (RECHOT  .LT.  51  THEN 
HOTCOL  =  1  ♦  (3iL0CATE(I,l))  +  RECHOT 
ELSE 

HOTCOL  =  3»iL0CATE(I,l)-l)  ♦  RECHOT-2 
END  IF 

C  DO  NOTION  CALCULATIONS 

IF  (RECHOT  .LT.  5)  THEN 
SYSHOT  =  HOTNS  (HOTCOL) 

RECNT  =  1 

40  IF  (RECDEG(RECNT,2)  .ST.  SYSHOT)  THEN 
IF  (RECNT  .EQ.  1)  THEN 
L0DE6  M.O 
LCHOT  »  0.0 

HIDES  =  RECDE6 (RECNT. 1) 

HIHOT  =  RECDES(RECNT,2) 

ELSE 

LODES  =  RECDES(RECNT-1,I1 
LOHOT  =  RECDE6(RECNT-1,2) 

HIDES  =  RECDES(RECNT,I) 

HIHOT  =  RECDES(RECNT,:) 

ENDIF 

VALUES(I)  =  LINTERP(SYSHOT. LOHOT, LOOES.HIHOT.HIDES) 

ELSE  IF  (RECNT  .EQ.  10)  THEN 
VALUES (I)  =  RECDEG (RECNT, 1) 

ELSE 

RECNT  *  RECNT  +  1 
GO  TO  40 
ENDIF 

C  DO  SPECTRUH  ANAL 

ELSE 
C 

DO  J  =  1,25 
RECNT  =  1 

DO  NHILE  (RECDEG(RECNT,2)  .LT.  HOTSPCIJ,!)) 

RECNT  =  RECNT  +  1 
END  CO 
C 
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IF  (RECNT  .EQ.  I)  THEN 
L0DE6  =  0.0 
LOHOT  =  0.0 
HICE8  =:  RECDE6(1,1) 

HIHOT  =  RECDEG(1,2) 

ELSE 

LODES  =  RECDES(RECNT-1,1) 

LOHOT  =  REC0E6(RECNT-1,2) 

HIDES  =  RECDES(RECNT,1) 

HIHOT  =  RECDE6(RECNT,2) 

END  IF 

IFER  =  LINTERP  (HOTSPC(J,I),LOHOT.i.uDEG.HinOT,HIDES) 
PRFSPC(J)  =  IFER  *  HOTSPC(',rtOTCOL) 

END  CO 
C 

IF  (PRFSPCIl)  .LI.  0.0)  THEN 
PRFSPCd)  =0.0 
END  IF 
PSHH  =  0.0 
DO  2060  J=2,25 
IF  (PRFSPC(J)  .LT.  0.0)  THEN 
PRFSPC(J)  =  0.0 
END  IF 

DIFF  =  HOTSPC(J,l)-HOTSPC(J-t,l) 

PSUH  =  PSUH  ♦  (PRFSPCIJ-1)+PRFSPC(J))»DIFF 
2060  CONTINUE 

VARIANCE  =  PSllH/2.0 

VALUES(l)  =  EFFEC (SORT (VARIANCE),  ERRLIH) 

END  IF 
END  IF 

r 

u 

2030  CONTINUE 
C 

50  TO  9099 

99?  STOP  ’  DATA  E.RRCfi  IN  SUBROUTINE  GETHOT’ 

9999  CONTINUE 
RETURN 
END 
C 
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C 

SUBROUTINE  OUTPUT (RONS, EFFECTIVENESS) 

C 

CHARACTERtlZ  NARES(S0,2) 

CHARACTERtl  TYPE (50, 2) 
lNTE6ERt2  LOCATE (50,2),  RONS, 1,J, LEAST 
REAL  VALUES (SO), EFFECTIVENESS 
C 

CONHON  /SYSTEN/  NAMES,  TYPE,  LOCATE,  VALUES 

r 

W 

NRITE  (*,120)  EFFECTIVENESS 
NRITE  (7,120)  EFFECTIVENESS 

120  FORMAT!///, 51, ’THE  EFFECTIVENESS  RELATIVE  TO  THIS  FEASIBLE  ’,/ 
k  ,51, ’DESIGN  CALM  NATER  PERFORMANCE  IS’,//,10I,F4.2) 

C 

LEAST  =  1 
C 

DO  1=1, RONS 

IF  ((TYPE(I,2)  .EQ.  ’E’l  .AND.  (VALUES(I)  .LT.  VALUES (LEAST))) 
k  THEN 

LEAST  =  I 
END  IF 
END  DO 
C 

NRITE  (*,130)  NAMES (LEAST, 2), VALUES (LEAST) 

NRITE  (7,130)  NAMES(L£AST,2),VALUES(LEAST) 

130  FORMAT!/, ST. ’EQUIPMENT/SUBSYSTEM  NITH  NORST  PERFORMANCE’, 

4  /,10T, ’NAME’, 27T,A12,/10T, ’RELATIVE  EFFECTIVENESS’. 5T,F4.2) 

C 

NRITE  (*,50) 

NRITE  (7,50) 

50  FORMAT  (/,’  SYSTEM  LIST’, /,10T,’ BLOCK’, 9X,’ ELEMENT’, 

4  5T , ' SER/PAR’ , 2T, ’ BLOCK/ELEMENT’ , 2T, ’ EFFECTIVENESS’ ) 

DO  I  =  1,S0NS 

NRITE(*.100)NANES(I,1),NANES(I.2),TYPE(I,1>,TYPE(I,2),VALUES(I) 
NRITE(7,100)NAHES(I,1),NAMES(I,2).TYPE!I,1),TYPE(I,2),VALUES(I) 
100  F0RMAT(10T,A12,2T,A12,3T,A1,11T,A1,8T,F4.2) 

END  DO 
RETURN 
END 
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C  ItMIttIttttttttttttttttIttttItItItlllltlIttHHtttItItItItllftItIt 

C 

REAL  FUNCTION  LINTERP  (TGTX.LOX.LOKY.HII.HnY) 

C 

REAL  T6T1(,L0X,L0XY,HIX,HIXY 
C 

LINTERP  =  LOXY+((HIXY-LOXY)/(HIJ-LOI))t(TeTX-LOXl 
C 

RETURN 

END 

C 

c 

c 

c 

REAL  FUNCTION  £FFEC(SI6HA,LIHIT) 

C 

REAL  CONST,  SIGNA,  INCREtt,  LINIT,  DELTA,  X 
C 

CONST  =  SORT (2f3. 14159) 

INCREN  =  LiniT/50.0 
X  =  -UNIT 

EFFEC  s  EXP(-(X»*2)/(2»(SIG«A»t2))>  /  (CONSTtSIGHA) 

r 

W 

DO  1*1,99 
X  *  X  +  INCREK 

EFFEC*EFFEC  ♦  2»£XP(-(X«2)/(2f(SIGnA»*2)))  /  (CONST»SISHA) 

END  DO 
C 

X  »  LIMIT 

EFFEC  =  EFFEC  +  EXP(-IXH2)/(2f (SIGBAf«2)))  /  (CONST*SI6«A) 

EFFEC  *  (EFFEC/2.0)  »  INCREN 

RETURN 

END 


